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1. Introduction 

A precise determination of the Kaon B-parameter is required to constrain the CKM unitarity 
triangle analysis. Bk is defined in terms of hadronic matrix elements which can be computed using 
lattice QCD. At present, the inclusion of dynamical quark effects is an essential requirement in 
these lattice calculations. Indeed, a remarkably good agreement has been found between several 
independent quenched determinations of Bk leaving the "quenching" effects as the largest un- 
certainty. A complete study of the systematic effects, other than quenching, has been performed in 
ref. [Q] and reviewed in this conference [Q]. There are indications that once dynamical quarks are 
taken into account, a major source of uncertainty on lattice results arises from the determination 
of the renormalization factors of the four-fermion operators |jl]|. In principle, this uncertainty can 
be completely eliminated by a non-perturbative renormalization procedure. Similarly, the renor- 
malization group (RG) running of the operator from hadronic scales up to high-energies, where 
perturbation theory can be safely applied, is best performed non-perturbatively. Here we report on 
the status of such a non-perturbative renormalization using Nf = 2 dynamical flavours. In particu- 
lar, we present preliminary results for the scale evolution of the four-fermion operator relevant for 
the determination of Bk in the context of tmQCD [ffl]. 

The theoretical description of the K° — K° oscillation is controlled, once high-energy scales are 
integrated through an operator product expansion procedure, by the matrix element (K Q \0 AS=2 \K°), 
where the four-fermion operator is defined as follows: 

AS = 2 = [jy„(l - y 5 )d) [sy M (l - y 5 )d] = <9 VV+AA - <9 VA+AV . (1.1) 

The strange and down quark fields are denoted by s and d respectively. The fi^-parameter is 
expressed in terms of the parity-even O vv +aa operator : 

_ (K°\O ww+AA \K ) 

B *= 8^2 • d- 2 ) 

3 r K m K 

Parity conservation ensures that the matrix element (K Q \O yA+A ^\K ), involving the parity-odd oper- 
ator, is identically zero. 

In lattice regularizations preserving chiral symmetry the operator V v+aa is multiplicatively 
renormalizable. This is not the case for Wilson fermions because in this case chiral symmetry is 
broken at non-zero lattice spacing. The renormalization of V v+aa is therefore more involved since 
mixing with four other dimension-6 operators has to be considered. On the other hand, discrete 
symmetries protect the parity-odd operator V a+av, so as to preserve multiplicative renormalization 
also in the case of Wilson fermions [js} §]. 

The inclusion of a "twisted mass" term in the fermionic action opens the way to improve 
the renormalization properties of Wilson fermions. The twisted-mass theory is related to standard 
QCD through an axial transformation of the quarks fields. By choosing appropriate formulations of 
tmQCD (see ref. [Q] for two of these formulations), it is possible to relate the QCD V v+aa operator 
to a partner V a+av in tmQCD where the property of multiplicative renormalization of this operator 
is preserved. This method allows the determination of Bk only via multiplicative renormalization, 
thus avoiding mixing with operators of wrong chirality. Another proposal in this direction uses 
chiral Ward-Takahashi identities to relate parity-odd and parity-even operators [j7|]. 
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Our strategy to renormalize the operator V a+av closely follows the one used in ref. ^ for the 
quenched case. The connection between the renormalization group invariant operator V a+av and 
its bare counterpart V a+av(§o) can be written in the following way: 



: lim ^VA+AV;s 

go^O 



O v 



r(x;go)- 



(1.3) 



The RGI operator is independent of the renormalization scheme and scale when the renormaliza- 
tion conditions are imposed at zero quark mass [g]. The renormalization factor ^v A +Av ;s (go) is 
scale-independent but depends on the scheme s (only through cutoff effects) and on the lattice 
regularization. These dependences are manifest when decomposing i2v A+ Av ;s (<?o) into: 



i^A+AV; S (§0) 



(1.4) 



The first factor on the r.h.s. of eq. ( |1.4[ ) controls the RG-running of the operator from the refer- 
ence scale /Vj n to an infinite scale. It is independent of the regularisation. The second factor, 
Z VA+AV . s (go,«|Umin)j relates the bare lattice operator to its continuum value at the hadronic scale 
jU m i n . This factor is therefore dependent on both the regularization and the scale. Both factors 



on the r.h.s. of eq. (1.4) depend on the renormalization scheme. In this report, we focus on the 



more expensive part of the renormalization procedure, that is, the computation of the contribution 
of the non-perturbative evolution function Z VA+ Av ;s (A t min) describing the running in the scale range 
1 GeV - 100 GeV. 



2. Renormalization group running of four-fermion operators 

Let us first define a renormalized four-fermion operator Or at a reference scale /I : 

R (x;n) = lim Z {go,aii) O{x;g ). 



(2.1) 



The running of the renormalized operator Or{}x) is controlled by its anomalous dimension Jo{g), 
defined as : 

M j-Or(xiii) = Yo(g) R (x;n) . (2.2) 

In mass-independent renormalization schemes [^] as those we consider here, the function Yo(g) 
only depends on the renormalized coupling g. The perturbative expansion of Jo{g) is given by 



ro(g) s ~°-g 2 (r$ ) + ri; ) g 2 + 7g ) g 4 + . 



(2.3) 



with y a universal coefficient. By combining eqs. (|2JJ) and ( |2.2| ), it is possible to relate the 
anomalous dimension of Or to its scale-dependent renormalization factor : 



7oOO)) =Jim UL^Z (go,an)\ Z o (g ,an) 1 . 



(2.4) 



The RGI operator is obtained upon the formal integration of eq. (2.2). It is given by 



6{x) = R (x;n) 



goo 

471 



i 0) /(2M 



exp • 



m ds (ro(g) rff r 



(2.5) 
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The integral in the r.h.s of eq. (2.5) describes the scale evolution of Yo(g)- The evolution function 



of the operator Or between the renormalization scale /I and an arbitrary scale \x' is given by: 

Uf(M) P{S) J a^o Z {go,aiJ.) 

The running of the renormalized four-fermion operator Or can therefore be performed by con- 
structing ratios of the renormalization factors Zq at different scales. 

Our renormalization schemes are defined in the Schrodinger functional (SF) formalism. This 
technique has been used to determine the scale evolution of physical quantities such as the strong 
coupling JTo| , 11] or the quark mass [12, 13]. These studies were performed both with and with- 
out dynamical quarks. In the case of four-fermion operators \^ the running was carried out in 
the quenched approximation. We regularize the theory on a lattice of physical size L 4 using stan- 
dard SF boundary conditions allowing to carry out simulations at zero quark masses. The SF is 
used as a mass-independent renormalization scheme; since the renormalization factors are flavour- 
independent, they can also be used to renormalize the B-parameters in the Kaon, D and B-meson 
sectors. The renormalization conditions are imposed at a scale pL equal to the IR cutoff 1 /L. 

Let us now concentrate on the case of the local parity-odd four-fermion operator: 

Ova+avW = ^ [(TOV^X^y^V^) + {ViYiiYs V2)(%Y» ¥4) + (¥2 <-» • (2-7) 

Four distinct valence flavours are used in the definition of the operator. The SF correlation functions 
used to extract the operator V a+avM are: 

% A ,r B ,r c ] (*d) = ^3 (^21 [T A ] ^ 45 [r B ] VA+AV (x) ^53 Fc]> , (2-8) 

where 6 and C' are the interpolating fields on the time boundaries (refer to [||] for a full explanation 
of the notations). Several choices of the Dirac matrices Fa,b,c are allowed. We will focus here on 
the particular choice Ta = Tb = Tc = 75. Note that a "spectator" valence quark 1//5 is used in 
eq. ( |2.8| ). It is useful to keep in mind that since the quarks are massless in this mass-independent 
renormalization scheme, flavour only enters through Wick contractions. 

The logarithmic divergences of the local operator V a+av (x) are isolated by dividing out from 
the correlator & the divergences coming from the boundaries and the external legs. This is obtained 
through the ratio: 

Kxo) = ^MgM, (2 . 9) 
/1 

where f\ is the boundary-to-boundary correlation function : f\ = — 1/ (2L 6 ) (^ 2 [ys] ^21 [Ys]) ■ The 
renormalized ratio Hr can be written as follows: 

h R (x ;n) =Z VA+AV (g ,an)h(x ;go) , (2.10) 

where the renormalization factor is fixed by imposing the renormalization condition : 

Z WA+AV (g ,an = l/L)h(x =L/2;g ) = h s (x =L/2;g ) , (2.11) 

so=o 
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i.e. at tree level Z VA+A v = 1- The renormalization condition is taken at the scale \i = \/L and 
therefore at fixed renormalized coupling u = g 2 (\/L). Note that as the quarks are massless, once 
the continuum limit a — > is taken, the only remaining scale is L. 

The running of the scale-dependent factor Z WA+AV (go,a/J.) is implemented in the SF formalism 
via the step scaling function (SSF), defined in the continuum as: 



f \ i- v / ir\ v i /r\ Z VA+AV (g(),fl:/2L) 

,(u) = lim£ VA+AV (M, a/L) , Z WA+AV (u,a/L) = — — - 

ti^o ^ A+AW {g ,a/L) 



. (2.12) 

m=0, g 2 (l/L)=n 



The SSF can be written in terms of the evolution function U : ct va+av (w) = Uy A+AW (l/2L, 1 /L) . The 
SSF is used to run the operator O va +av between two scales differing by a factor of two. By iterating 
this procedure the running of the operator can be performed over a large range of scales. 

3. Non-perturbative study of the step scaling function 

The computation of the SSF is performed with Nf = 2 flavours of 6(a) improved Wilson 
fermions. We evaluate L WA+AV (u,a/L) at six values of the renormalized coupling u (labelled, in 
increasing order, u\,...,u6). At each of these couplings, we consider three lattice resolutions L/a = 
6, 8, 12 to extrapolate our data to the continuum limit. The unquenched configurations employed in 
our computation have been previously used in the study of the quark mass renormalization flT3[ ] (the 
description of our simulation setup can be found in this reference). In fig.|l| we present the status 
of the ongoing determination of L WA+AV (u,a/L). Data for some of the simulation points has not yet 
been included as it is still being generated. Moreover, the statistical errors of our preliminary data 
in ng.[l] will reduce when the complete set of configurations will be considered. The integrated 
autocorrelation times are included in the error estimate. We observe that the autocorrelations grow 
when increasing the coupling and when approaching the continuum limit. 

The 0(a) improvement of the dimension-six operator V a+av has not been implemented. Al- 
though a Symanzik improvement program is possible, the mixing of V a+av with several dimension- 
seven operators makes it unpractical. For each of the couplings u, the continuum limit of Tva+av (u,a/L) 
should therefore be taken through a linear extrapolation. This is illustrated in fig.[[] in those cases 
where three L/a resolutions are already available. As our data seems to show rather small cutoff 
effects, we have also tried to fit L WA+A ^(u,a/L) to a constant (we perform a weighted average and, 
somehow abusively, we refer to it as a "constant fit"). In the case of the couplings ul,u3,u5 and 
w6 this fit was performed by discarding the L/a = 6 data which, being far from the continuum, is 
expected to contain large cutoff effects. By comparing the linear and the constant fit, we observe 
good agreement of the extrapolated values. A more refined analysis of the continuum extrapolation 
will be undertaken once our complete set of data will be available. 

In a tentative study of the quality of our data, we present in fig. || the SSF <7 V a+av(w)- The Nf = 2 
data is compared to the quenched one from ref. [||]. The same renormalization scheme and fitting 
procedures is used in both the Nf = 2 and Nf = data (in particular, also in the quenched case we 
perform a constant fit to the continuum). 1 In the strong coupling regime, u ~ 3.5, we observe 

1 Note that this comparison is only intended to study the quality of our data. It is indeed improper to compare Nf = 2 
and quenched physical results at this stage since the respective renormalized couplings are not taken at the same physical 
scale. 
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Figure 1: Continuum extrapolation of the SSF E V a+av(m,<V^) at fixed renormalized couplings m. The empty 
circles correspond to the extrapolated value obtained through a linear fit and the diamond to the one obtained 
via a fit to a constant. Results are preliminary. 



a similar pattern in the Nf = 2 and Nf = data when comparing the value of (J VA+AW (u) obtained 
through a linear and a constant extrapolation to the continuum: in both Nf = 0,2 cases the linear 
extrapolation points lie above the constant-fit ones. Due to large statistical errors, the Nf = 2 data 
shows a better agreement between linear and constant fit. In fig.||, we also plot the perturbative 
expressions of a VA+ Av(w) in both the Nf = 2 and iVf = cases. These expressions were computed 
at next-to-leading order (NLO) in ref. [14]. In our chosen renormalization scheme, we observe a 
fairly good agreement between the perturbative curve and the non-perturbative data in the small 
coupling region and some signs of deviations in the strong coupling regime. We have considered 
nine different renormalization schemes (for the definitions of these schemes, refer to [|8, 14]). The 
empirical criterion to identify the more appropriate schemes is to consider those having a small 
NLO term, of the same sign as the LO one, in the perturbative expansion of cj va+ av(w)- We have 
checked that our best available scheme is indeed the one of eq. (|2T9|). 



Conclusions 

We have presented the status of the computation of the non-perturbative RG running of the 
four-fermion operator V a+av using lattice QCD with two dynamical quarks. This calculation will 
soon allow us to determine the universal renormalization factor Z V a+av;s 

0) in eq. Q. The second 

factor on the r.h.s of this equation is simpler to determine, compared to Zv A+AV . 6 (jU), as it depends 
only on a single scale ll. This determination will allow to complete the renormalization of the 
operator V a+av- As our renormalization scheme is flavour-independent, the same renormalization 
factors can be used to determine the B-parameters in the strange, charm and beauty sectors. 
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Figure 2: The step scaling function OVa+av( m ) (discrete points) is compared to the LO and NLO perturbative 
results. In order to evaluate the quality of our data we compare our Nf — 2 results with quenched data from 
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